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Abstract 

We explore two possible ways to enhance the critical temperature Tc in the dilute magnetic 
semiconductor Mno.08Gao.92As. Within the context of the double-exchange and RKKY pictures, 
the ternary alloys Mn2;Cro.o8-a:Gao.92As and Mno.o8Gao.92Asj/Pi_y might be expected to have 
higher than the pseudobinary Mno.08Gao.92As. To test whether the expectations from model pic- 
tures are confirmed, we employ linear response theory within the local-density approximation to 
search for theoretically higher critical temperatures in these ternary alloys. Our results show that 
neither co-doping Mn with Cr, nor alloying As with P improves Tc- Alloying with Cr is found to 
be deleterious to the Tc. Mno.o8Gao.92Asj;Pi_j, shows almost linear dependence of Tc on y. 

PACS numbers: 75.50.Pp, 75.30.Et, 71.15.Mb 
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Searching for spintronic materials with high Curie temperature Tc has attracted a lot 
of interest recentl}*ii^ because of the promising future of these materials. In spintronics 
technology, by including both magnetic and electronic properties in the devices, we can 
manipulate the electronic spin to the same degree as what we do to the electronic charge 
in conventional semiconductor devices. It can greatly improve the performance of current 
digital information processing and storage technology. Currently, one of the central issues for 
spintronics application is the development of carrier-mediated ferromagnetic semiconductors 
that are magnetic above room temperature. Dilute magnetic semiconductors (DMS), i.e. 
semiconductors doped with low concentrations of magnetic impurities (such as Mn, Co, Cr), 
are generally thought to be good candidates for these spintronics applications. Gai_a;Mna;As 
is one of the most widely studied DMS, because it is one of the few materials where it is 
generally agreed that the magnetism is carrier-mediated. 

Currently, much effort has been expended to increase in Gai_a;Mna;As and related DMS 
materials. Tc has risen steadily in Gai-^jMn^As, to Tc ~ 180K in very thin films, carefully 
annealed to eliminate Mn interstitials without at the same time allowing MnAs precipitates 
to fornt^^. However, there is still a large gap between the highest Tc obtained to date and 
a technologically useful working temperature. 

A systematic understanding of the ferromagnetism in DMS systems is necessary for the 
design and optimization of the materials. However, the theoretical understanding of this 
system is still not very clear, although several models were proposed. One generally used 
model is based on the assumption that Mn local magnetic moments interact with each other 
via RKKY-type interactions, the Mn at the same time providing a source of free holea^. 
This picture can explain the temperature dependence of the magnetization and the hole 
resistivity of the DMS, although there is a growing consensus that many of the claims of that 
paper were artifacts of the assumptions in the model. On the other hand, Akai^ first used 
the local spin-density approximation (LSDA) to estimate Tc within the Coherent Potential 
Approximation (CPA) in (In,Mn)As; he argued that a double exchange mechanism was a 
more appropriate description of the magnetism than the pd exchange assumed by Dietl. 

In addition, there is a fair amount of research that indicates the ferromagnetism in DMS 
systems are strongly affected by the inherent disorder. In our recent work^ we included 
various effects of disorder, and found that there is a rather strict upper limit of around 250K 
for Tc in disordered Mn:GaAs systems (as computed within the LDA), which implies a rather 
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unpromising future for the simple Mn:GaAs alloy. So, it is becoming an important issue to 
find possible ways to optimize the magnetic properties of DMS. The possible means may 
include co-doping the magnetic elements, change of host materials and magnetic elements, 
or spacial confinement of the magnetic dopants, etc. We use the Zener double-exchange 
picture for guidance, for reasons that will be clear below. For detailed analysis we employ 
to the same LDA linear response techniques as in Ref.— . 



FIG. 1: Cartoon of the Zener model coupling levels on two neighboring sites. Top illustrates 
the majority and minority levels on each site together with the model parameters (intra-atomic 
exchange splitting J, effective interatomic hopping matrix element V , and canting angle 9). Bottom 
left depicts the four levels when spins are aligned parallel; bottom right depicts the levels when 
spins are aligned antiparallel. 

Here we consider two strategies to make alloys related to Mn^Gai_^As, that can have 
higher Tc. We based both strategies on the double exchange /super exchange picture^ depicted 
in Fig. H The four levels are described by the hamiltonian 
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with U the spinor rotation matrix for the second site: 
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The eigenvalues are given by the quartic equation^ 

e2 = J"^ + ±2VJcos{e/2). 



(3) 



When only the lowest state is filled, the parallel alignment dominates (Zener double ex- 
change). The most salient other characteristic of the Zener double exchange energy is its 
unusual non-Heisenberg dependence on angle, varying as |cos(0/2)|. When the first two 
states are filled the energy vanishes for parallel spin alignment, and the the system is stabi- 
lized in the AFM alignment (Anderson superexchange), the energy varying as —cos{9). As 
a function of filling, there is a transition from FM stabilization to AFM. The maximum FM 
stabilization occurs at half-filling. 

In the Mn^^Gai-xAs and Crx^Gai^^As alloys, the exchange is mediated by the partially 
filled TM-derived t2 impurity band. In the dilute limit, the Mn-derived ^2 level has threefold 
degeneracy and falls at VBM+0.1 eV, while the Cr-derived level falls at VBM+0.3 eV. These 
levels broaden into an impurity band as x increases. In reality this band is composed of an 
admixture of TM d states and effective-mass like states of the hosl^^^, mostly derived from 
the top of the host valence band. 

To the extent that the model correctly describes exchange in Mn^^Gai-^^As, it is possible 
to optimize T^. of a specific DMS system by varying the carrier concentration to set this 
impurity band as close to half filling as possible. One possible way to do this is to make a 
pseudoternary alloy Mn2;CryGai_i;_j/As, where the (Mn,Cr,Ga) atoms share a common (fee) 
sublattice. Cr has one electron to fill the 3-fold majority t2 level to 1/3 filling. Mn has two 
electrons, so the majority t2 level should be 2/3 full. By alloying Mn and Cr, we might 
obtain a compromise 1/2 filling that will increase Tc. 

Another possibility is to vary the anion in the III-V series. From model viewpoints (and 
confirmed by the LDA), there are two competing effects from the anion. The important 
quantity is the position of the transition-metal d level relative to the top of the valence 
band (determined largely by the anion p level). When the anion is heavy (e.g. Sb) the p 
is shallow; the Mn d sits mostly below the valence band, and the ^2 derived impurity band 
is mostly composed of the effective-mass like host states; thus the exchange is effectively 
mediated the holes near the valence band top, as the RKKY picture assumes. If the Mn d 
is very deep, the coupling between the Mn and holes also weak, so the interatomic exchange 
is small. As the anion becomes lighter (Sb^As^P^N) the anion p deepens, so that finally 
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the Mn d falls above the valence band and forms a deep level. In this limit the coupling 
should be well described by double-exchange. For Mn in GaN, the t2 level is approximately 
at midgap. It couples strongly to nearest neighbors, but is short ranged. At least in the 
dilute case, the near-neighbor exchange interactions are mostly ferromagnetic and rather 
strong, but they decay rapidly with distance^. But also in the dilute case, there are few 
nearest neighbors and more distant interactions are necessary for long-range order to exist, 
as without them there is no percolation path. The net result is that Tc is predicted to be 
quite low in Mn^Gai_j;Ni2i. The anion with the highest T^. will adopt some compromise 
between a very deep valence band (GaN) and a very shallow one (GaSb). With this picture, 
and the expectation that the Mn d level in MnGaP is perhaps already too high, we might 
be able to make a ternary MnGaPa;Asi_^ where x can be adjusted to tune the position of 
the valence band, and thus optimize Tc. 

First principles calculations provide a way for evaluation of the exchange interactions 
and Tc in a DMS system of a specific geometry. We adopt here the same scheme as pre- 
sented in Ref.^, to study the ferromagnetic trends for the ternaries Mna;Cro.o8-xGao.92As and 
Mno.osGaPyAsi-y. Self-consistent calculations of 200-atom special quasirandom structures'^, 
or SQS, were performed using the method of Linear-Mufhn-Tin Orbitals (LMTO) within 
the local spin density approximation (LSDA) and the Atomic Spheres Approximationi^ii^. 
The exchange interactions were computed within the long-wave approximation using a linear 
response Green's function technique, which maps the total energy onto a Heisenberg form. 
The interatomic Heisenberg exchange parameters are determined from the linear response 
expression 

4 = ^ r ^^Tr, {P, . {tI -Tl^ + ry. Tl) . P, } (4) 

where i and j are indexes of atoms with the on-site perturbation having a form P = Pz = 
(Tf ^ - Tf ^)/2 = (PT - pi)/2, as described in Ref.i^. P is a diagonal matrix whose rows 
and columns correspond to the orbitals at a particular site {spd in the present case); P 
is closely related to the majority- minority spin-splitting. tJj ■ Tj^ essentially the transverse 
susceptibility. Magnetization at finite temperature and Tc were calculated using the Cluster 
Variation Method (CVM) of KikuchiiSi, generalized to the continuous rotational degrees of 
freedom of the (classical) Heisenberg model, which predicts Tc to within about 5% of results 
calculated by full spin-dynamics simulations in these DMS alloys^. For both Gao,92Mno,o8As 
and Gao.92Cro.08As the CVM predicts Tc of about 280K. We will show that, although these 
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two systems are intriguing based on the simple pictures we carry around about them, neither 
alloy is predicted to have magnetic properties more favorable than Mn^^Gai-ajAs at a; = 0.08. 

Before addressing T^. in the pseudoternary Mn^CryGai_^_j^As, we consider a simple model 
calculation to see to what extent the LDA exchange interactions approximately parallel the 
double-exchange picture. Starting from Mno.08Gao.92As SQS configuration, we artificially 
shift the Fermi level of this system to mimic additional doping. In this model dopant 
electrons are assumed to have no effect but change the filling of the valence band, which 
enables us to isolate this effect from other complexities that occur in the real system, but 
which are missing in the double-exchange model. 
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FIG. 2: Dependence of Tc Fermi level shift and the corresponding change in band occupancy. 
Calculated MFA and CVM Tc are presented as a function of Fermi level shift and deviation from 
charge neutrality. 

Fig. 121 shows the change in Tc due to the Fermi level shift AEp (and corresponding hole 
doping), in a 200-atom SQS structure Ga92Mn8Asioo- As predicted by both the RKKY 
picture and the double-exchange picture, increasing Ep reduces the number of holes and 
reduces Tc. Conversely, with the initial introduction of more holes [AEp < 0) increases 
as both RKKY picture and double-exchange models predict. However, fig. El shows that Tc 
reaches a maximum for AEp ~ —0.20 eV; further shifts reduce Tc. Observing the Fermi 
level shifts in the DOS we can confirm that position of Ep at optimum Tc corresponds to 
Ep falling in the middle of the Mn ^2 band, indicating half filling of the t2 level. This can 
be determined independently by counting the number of holes added as a function of AEp. 
AEp ~ —0.20 eV corresponds to the addition of approximately 4 holes in this system of 
8 Mn atoms, or 0.5 hole/atom. The ^2 level with initially 2 electrons {AEp = 0) changes 
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to 1.5 electrons near AEp = —0.2, which is just half of the number needed to fill the 
level. This confirms that the double-exchange picture more closely reflects the ab initio 
LDA calculations than does RKKY, as Akai first proposed^. 

One possible way to introduce holes in a real system is to add another dopant, e.g. Be, as 
has been done in (5-doped MnGaAs^^ However, Be doping appears to catalyze the formation 
of other deleterious defects (probably Mn interstitials, which being donors, bind to Be). 
Another way is to alloy Cr with Mn. In the simplest view we can mimic this effect by a 
"virtual crystal" where we use a hypothetical atom X with atomic number 24.5. The 
evaluated for this hypothetical Ga92X8Asioo SQS structure is dramatically larger than Tc 
for Ga92Mn8Asioo or Ga92Cr8Asioo, as we found before. For this 8% doping case, the T^. is 
calculated to be 710K with the MFA and 390K within the more reliable CVM. 




FIG. 3: Dependence of Tc on y in Mn^CryGai^x-y-^s. Tc is computed both within the MFA and 
the CVM. The noise in the data occurs because Tc depends on the specific choice of quasi random 
configuration. For the same y, there can be a deviation between different SQS structures as much 
as 50K. 



Finally, we consider a physically realizable configuration, Mn^CiyGai-x-y-^^- We now 
construct a series of three-body 200-atom SQS structures, varying x for fixed x + y=8%. 
Tc was calculated for all a variety of configurations with different x. The calculated Curie 
temperature depending on x are shown in Fig. El The picture that emerges is, unfortunately 
very different from the virtual crystal. Co-doping actually decreases T^. We can see although 
CVM and MFA Tc curves have a little different shape, they both predict Tc to decrease near 
X = 1/2, in contradistinction to the virtual crystal. The reason can be traced to the fact that 
Cr and Mn have their impurity levels at different energies. Consider the limit of infinitely 
narrow impurity bands: there is a Cr-derived level at VBM-l-0.3 eV and a Mn-derived level 



at VBM+0.1 eV. The single electron in the Cr-derived level falls to the Mn level, thus 
emptying out the Cr t2 state and completely filling the Mn t2 state. In the double-exchange 
picture, the system will be antiferromagnetic. In the alloy, these levels broaden into bands. 
However the Mn and Cr like states derive from sufficiently different energies, so that rather 
than forming a single wide impurity band made of equal admixtures of Cr and Mn, the band 
is composed of a lower part more of Mn character and an upper part dominated by Cr. If 
the bandwidth is not sufficiently large to counter this effect, the net result is that Tc{x) in 
the true ternary apparently bows oppositely to what occurs in the idealized case. 



FIG. 4: Energy resolved Partial DOS of Mn and Cr ions in Ga92Mn4Cr4Asioo, showing similar 
filling as single element doping case. 

This is reflected in the partial DOS of TM in (Mn,Cr):GaAs(Fig. Ej). The DOS of the 
for each individual TM ion reflects the DOS in the respective pseudobinary alloy: the Mn 
t2 near Ep is centered slightly below the Cr t2, showing the localized behavior of magnetic 
dopants. Some broadening is evident (particularly in the tail of the Mn t2 above Ep, but 
owing to the dilute alloy concentration the coupling between t2 levels is not sufficient to 
broaden it into a homogeneous band. The lack of this distinction makes the virtual crystal 
approximation an unsuitable description of the real alloy. When both levels are present 
the deeper Mn t2 (2/3 full in Mn:GaAs) becomes fuller still, while the Cr t2 (1/3 full in 
Cr:GaAs) becomes more empty. Hence, the co-doping actually drive the filling opposite to 
the desired direction. The band is actually adversely affected, resulting in a decrease in 
by co-doping. 
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It is notable that an RKKY-like picture^, which describes exchange interactions purely in 
terms of spin-split host bands whose splitting is determined by the local magnetic moments, 
will also miss this effect. This establishes in another way the fact that exchange interactions 
within the LDA more closely reflect the Zener double-exchange model than an RKKY-like 
model. 

Our other strategy to increase Tc is to consider the ternary Mno.o8Gao.92AsyPi_j/, as 
we described above. In this calculation, SQS configurations are also generated through 
independent sublattices of cations (Mn,Ga) and anions (As,P). the Mn concentration is 
fixed to be 0.08% as before. Tc for different P concentrations are studied. For the lattice 
constant, we use virtual crystal approximation with modification based on Ref.— . The 
figure shows that the net effect (deviation from linear behavior) is small, but in any case 
Tc bows downward with y, with Tc suppressed by ~50K for y = 0.5. This reduction can be 
understood as a suppression, through disorder, of the number of paths that connect different 
magnetic sites, and thus the magnetic exchange between them. 

600 1 ■ 1 ■ 1 ■ < ■ 1 
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FIG. 5: Curie Temperature for Mno.o8Gao.92AsyP function of y 

In conclusion, based on the LSDA, we found out the Curie Temperature of DMS can 
be improved by modifying the band filling, as both double exchange and RKKY models 
predict. We show that the LSDA results more closely parallel the double exchange than the 
RKKY picture, and predict that optimal Tc occurs when the impurity ^2 band is half full. 
However, when we attempt to modify the concentration in a real crystal like Mn and Cr 
co-doping GaAs, Tc is actually suppressed. For a successful application of this rule, it will 
be necessary to find another way to to tailor the band occupancy in order to achieve half 
filling. We also found that admixing of As and P on the anion sublattice does not improve 
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